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FAILURE CRITERIA FOR FRP LAMINATES 
IN PLANE STRESS 


Carlos G. Davila and Pedro P. Camanho 


Abstract 

A new set of six failure criteria for fiber reinforced polymer laminates is 
described. Derived from Dvorak’s fracture mechanics analyses of cracked 
plies and from Puck’s action plane concept, the physically-based criteria, 
denoted LaRC03, predict matrix and fiber failure accurately without 
requiring curve-fitting parameters. For matrix failure under transverse 
compression, the fracture plane is calculated by maximizing the Mohr- 
Coulomb effective stresses. A criterion for fiber kinking is obtained by 
calculating the fiber misalignment under load, and applying the matrix 
failure criterion in the coordinate frame of the misalignment. Fracture 
mechanics models of matrix cracks are used to develop a criterion for 
matrix in tension and to calculate the associated in-situ strengths. The 
LaRC03 criteria are applied to a few examples to predict failure load 
envelopes and to predict the failure mode for each region of the envelope. 
The analysis results are compared to the predictions using other available 
failure criteria and with experimental results. Predictions obtained with 
LaRC03 correlate well with the experimental results. 


Introduction 

The aim of damage mechanics, the mathematical science dealing with quantitative descriptions of the 
physical events that alter a material when it is subjected to loads, is to develop a framework that de- 
scribes the material response caused by the evolving damage state. The greatest difficulty in the develop- 
ment of an accurate and computationally efficient numerical procedure to predict damage growth has to 
do with how to analyze the material micro-structural changes and how to relate those changes to the 
material response. Several theories have been proposed for predicting failure of composites. Although 
significant progress has been made in this area, there is currently no single theory that accurately predicts 
failure at all levels of analysis, for all loading conditions, and for all types of fiber reinforced polymer 
(FRP) laminates. While some failure theories have a physical basis, most theories represent attempts to 
provide mathematical expressions that give a best fit of the available experimental data in a form that is 
practical from a designer’s point of view. 

To the structural engineer, failure criteria must be applicable at the level of the lamina, the laminate, 
and the structural component. Failure at these levels is often the consequence of an accumulation of 
micro-level failure events. Therefore, an understanding of micro-level failure mechanisms is necessary in 
order to develop accurate failure theories. 



The recent World Wide Failure Exercise 1 ' 7 (WWFE) conceived and conducted by Hinton and Soden 
provides a good assessment of the status of currently available theoretical methods for predicting 
material failure in fiber reinforced polymer composites. The recently published comparison of the 
predictions by the WWFE participants with experimental results indicates that even when analyzing 
simple laminates that have been studied extensively over the past 40 years, the predictions of most 
theories differ significantly from the experimental observations 3 . 

The uncertainty in the prediction of initiation and progression of damage in composites has led to the 
undertaking of an effort at the NASA Langley Research Center to revisit existing failure theories, assess 
their capabilities, and to develop new theories where necessary. The present paper describes a newly 
developed set of six nonempirical criteria for predicting failure of unidirectional FRP laminates. All of 
the calculations presented are at the lamina level with plane stress assumptions. First, the motivation for 
non-empirical failure criteria for matrix shear and transverse tension and compression is established by 
examining several existing criteria. Next, new matrix failure criteria for matrix tension and compression 
are proposed. Then, a new fiber kinking failure criterion for fiber compression is developed by applying 
the matrix failure criteria to the configuration of the kink. The resulting set of six failure criteria is 
denoted LaRC03 and are summarized in the Appendix. Finally, some examples of failure envelopes are 
presented. The predicted results are compared to the results of other available failure criteria and with 
experimental results. 

Strength-Based Failure Criteria 

Strength-based failure criteria are commonly used with the finite element method to predict failure 
events in composite structures. A large number of continuum-based criteria have been derived to relate 
internal stresses and experimental measures of material strength to the onset of failure. One detailed 
discussion on the failure criteria proposed by Hashin, Puck, Sun and several others has been presented by 
Paris. 8 Paris also discusses the ad hoc nature of the formulation of most strength-based criteria. In the 
following sections, the Hashin criteria 9, 10 are briefly reviewed, and improvements proposed by Sun 11 and 
Puck 5,6 over Hashin’ s theories are examined. 

Hashin Criteria 2D (1980) 

Hashin established the need for failure criteria that are based on failure mechanisms. In his 1973 
proposal, he used his experimental observations of failure of tensile specimens to propose two different 
failure criteria, one related to fiber failure and the other related to matrix failure. The criteria assume a 
quadratic interaction between the tractions acting on the plane of failure. In 1980, he introduced fiber and 
matrix failure criteria that distinguish between tension and compression failure. Given the difficulty in 
obtaining the plane of fracture for the matrix compression mode, Hashin used a quadratic interaction 
between stress invariants. Such derivation was based on logical reasoning rather than micromechanics. 
Although the Hashin criteria were developed for unidirectional laminates, they have also been applied 
successfully to progressive failure analyses of laminates by using in-situ unidirectional strengths 12 to 
account for the constraining interactions between the plies. The two-dimensional versions of the failure 
criteria proposed by Hashin in 1973 and 1980 are summarized in Table 1. 

Numerous studies conducted over the past decade indicate that the stress interactions proposed by 
Hashin do not always fit the experimental results, especially in the case of matrix or fiber compression. It 
is well known, for instance, that moderate transverse compression (< 7 22 < 0) increase the apparent shear 
strength of a ply, which is not well predicted by Hashin’s criterion. In addition, Hashin’s fiber 
compression criterion does not account for the effects of in-plane shear, which reduce significantly the 
effective compressive strength of a ply. Several researchers have proposed modifications to Hashin’s 
criteria to improve their predictive capabilities. Modifications proposed by Sun and Puck are described 
below. 
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Improved Criteria for Matrix Compression 

Sun et al. 11 proposed an empirical modification to Hashin’s 1973 criterion for matrix compression 
failure to take into account the beneficial role that compressive a 22 has on matrix shear strength. The 
criterion is: 
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where 7j is an experimentally determined constant and may be regarded as an internal material friction 
parameter. The denominator S L -7ja 2 2 can be considered an effective longitudinal shear strength that in- 
creases with transverse compression a 22 (see sign convention in Fig. 1). As in Hashin’s theories, the 
angle of the fracture plane is not calculated. 



Figure 1. Fracture of a unidirectional lamina subjected to transverse compression and in-plane shear. 


Table 1. Hashin Criteria 9, 10 for Plane Stress 


Matrix Failure Matrix tension, <j 22 >0 



Fiber Failure Fiber tension, <j u > 0 



Matrix compression, g 2 2 < 0 
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Fiber compression, Gn < 0 



where, 

FI = Failure Index, where the subscripts M and F indicate matrix and fiber failure, respectively 
Oil = components of normal and shear stresses, with i,j =1,2 

X T , X° = strength in fiber direction under tension and compression, respectively (Longitudinal strength) 

Y r , Y 0 = strength normal to the fiber direction under tension and compression, respectively (Transverse strength) 
S L and S T are the longitudial and transverse shear strengths, respectively. 


Puck’s Action Plane proposal 5 represents the beneficial influence of transverse compression on matrix 
shear strength by increasing the shear strength by a term that is proportional to the normal stress <j n acting 
at the fracture plane shown in Fig. 1. In this formulation, the matrix failure criterion under transverse 
compression is: 
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where T t and T l are the shear stresses acting on the fracture plane defined in Fig. 1 . The direct 
contribution of a 2 2 has been eliminated by assuming that the initiation of fracture is independent of the 
transverse compressive strength. Internal material friction is characterized by the coefficients rf and rf \ 
which are determined experimentally. 

A key element to Puck’s proposal is the calculation of the angle of the fracture plane, a ; shown in Fig. 
1 . Puck determined that matrix failures dominated by in-plane shear occur in a plane that is normal to the 
ply and parallel to the fibers (oc=0). For increasing amounts of transverse compression, the angle of the 
fracture plane a changes to about 40°, and increases with compression to 53°±2° for pure transverse 
compression. In the WWFE, Puck’s predicted failure envelopes correlated very well with the test results . 6 
However, Puck’s phenomenological approach uses several material parameters that are not physical and 
may be difficult to quantify without considerable experience with a particular material system. 


LaRC03 Criteria for Matrix Failure 

In this section, a new set of criteria is proposed for matrix fracture that is based on the concepts 
proposed by Hashin and the fracture plane concept proposed by Puck. In the case of matrix tension, the 
fracture planes are normal to the plane of the plies and parallel to the fiber direction. For matrix 
compression, the plane of fracture may not be normal to the ply, and Hashin was not able to calculate the 
angle of the plane of fracture. In the present proposal, Mohr-Coulomb effective stresses 13 are used to 
calculate the angle of the fracture plane. 

Criterion for Matrix Failure Under Transverse Compression (q^ O) 

The Mohr-Coulomb (M-C) criterion 13 is commonly used in applications where fracture under tension 
loading is different from fracture under compression loading, such as in soil mechanics or in the fracture 
of cast iron. The application of the M-C criterion to multiaxial failure of epoxy resins was studied by 
Kawabata 14 based on correlation with his own test results. While studying the failure of chopped glass- 
fiber/epoxy mat laminates under confining pressures, Boehler 15 found the Tsai-Wu criterion to be 
inadequate, and formulated a shearing criterion based on the M-C criterion that fit his experimental 
measurements well. Taliercio 16 used the M-C criterion within a nonlinear micromechanical model to 
predict the macroscopic strength properties of fiber composites. 

The M-C criterion is represented geometrically by the diagram illustrated in Fig. 2. The Mohr’s circle 
represents a state of uniaxial compression. The angle of the plane of fracture a 0 is set in this example at 
53°, which is a typical fracture angle for composites under transverse compression loading (Puck 5 ). The 
line AB is the tangent to Mohr’s circle at A and is called the Coulomb fracture line. The M-C criterion 
postulates that in a state of biaxial normal stress, fracture occurs for any Mohr’s circle that is tangent to 

the Coulomb fracture line. The effective stress r e f f is related to the stresses z T and <J n acting on the 

fracture plane by the expression T e j r = z T + r/<J n . In the literature, tan -1 ( 77 ) is called the angle of internal 

friction and it is assumed to be a material constant. When rj= 0, the M-C criterion is equivalent to the 
Tresca condition 13 . 

DiLandro 17 explains the role of internal friction on the strength of carbon-fiber composites by noting 
the absence of chemical bonds between fiber and matrix, and that adhesion is attributed to Van der 
Waal’s interactions. When subjected to an external load, the shear slipping of the two phases is prevented 
until the shear stress at the fiber-matrix interface reaches a limiting value. DiLandro also notes that 77 is 
an empirical factor that encompasses all chemical-physical interactions, including the thermal residual 
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shrinkage of the matrix around the fiber. Larson 18 examined the relative effects of interfacial friction and 
roughness on the length of interfacial sliding which proceeds from the tip of an impinging fracture 
oriented perpendicular to the interface. According to Larson, sliding is key to the cracking behavior of 
fibrous brittle matrix composites in that it affects the stress concentration on the fibers, the matrix crack 
spacing, and, therefore, the global toughness of a composite material. 



-1 

tan(2a 0 ) 


Figure 2. Mohr’s circle for uniaxial compression and the effective transverse shear. 


In general, the fracture plane can be subjected to transverse as well as in-plane stresses, in which case 
the effective stresses must be defined in both orthogonal directions as shown in Eq. 3. 



+ n T cr n ) 

+ 7 l L(T n) 


(3) 


where the terms if and rj are referred to as coefficients of transverse and longitudinal influence, 
respectively, and the operand (x) = x if x > 0; otherwise (x) = 0. Matrix failure under compression 

loading is assumed to result from a quadratic interaction between the effective shear stresses acting on 
the fracture plane. The failure index for a failure mode is written as an equality stating that stress states 
that violate the inequality are not physically admissible. The matrix failure index ( FI m) is 


LaRC03 #1 


fi m = 


f t t 

L eff 
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where S T and S ! Ls are the transverse and longitudinal shear strengths, respectively. The subscript M 
indicates matrix failure. The subscript is indicates that for general laminates, the in-situ longitudinal 
shear strength rather than the strength of a unidirectional laminate should be used. The constraining 
effect of adjacent plies substantially increases the effective strength of a ply. It is assumed here that the 
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transverse shear strength S T is not subjected to in-situ effects, which is an assumption that deserves to be 
studied further. The calculation of in-situ strength is discussed in a later section of this paper. 

The stress components acting on the fracture plane can be expressed in terms of the in-plane stresses 
and the angle of the fracture plane, #(see Fig. 1). 

<j n = cos 2 a a 22 

z T = -sin# cos# <j 22 (5) 

t l = t u cos # 

Using Eqs. 3 and 5, the effective stresses for an angle of the fracture plane ^between 0° and 90° are 

z T eff = (- a 22 cos # (sin a-ij T cos #)) 

\ / ( 6 ) 

T eff = ( C ° S a (| T 12 | + r l L °22 C0S a )) 


Calculation of Coefficients rf and if and Strength S T 

The coefficients of influence rf and rf are obtained from the case of uniaxial transverse compression 
( c 22 < 0, t, 2 = 0 ). At failure, the in-plane compressive stress is equal to the matrix compressive 

strength, o 22 = -Y C - The effective transverse shear stress at failure is 

z T eff = S T - Y c cos a ( sin a - rf cos a) (7) 


Under uniaxial transverse compression, fracture occurs at a fracture angle a 0 that maximizes the effec- 
tive transverse shear. Taking the derivative of the transverse shear stress, Eq. 7, with respect to a gives 


d z 


eff 


da 


= 0 = - sin a 0 ( sin a 0 - rf cos a 0 ) + cos a 0 ( cos a 0 +rj T sin a 0 ) 
=^> cos 2 a 0 - sin 2 a 0 +2 rf sina 0 cosa 0 = 0 


( 8 ) 


Solving Eq. 8 for rf gives 
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-1 


tan 2a n 


( 9 ) 


Puck 5 determined that when loaded in transverse compression, most unidirectional graphite/epoxy 
composites fail by transverse shear along a fracture plane oriented at a 0 = 53°±2°. Therefore, the 
coefficient of transverse influence is in the range 0.2\<rf<036. Note that if the fracture plane were 
oriented at #^=45°, the coefficient of transverse influence would be equal to zero. 

The transverse shear strength S T is a material property that is difficult to measure experimentally. 
However, substituting Eq. 9 into Eq. 7 gives an expression relating the transverse shear strength to the 
transverse compressive strength: 


S T = Y c cos a 0 ( sina 0 + 


cos a 


tan 2 a, 




( 10 ) 
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For a typical fracture angle of a 0 =53° gives S T =0.378 as was shown in Fig. 2. Note that in some 
textbooks, the transverse shear strength is often approximated as S T =0.5 which implies from Eqs. 7 
and 9 that the fracture plane is at a 0 =45° and that rf=0. Also note that with this approximation, Hashin’s 
1980 two-dimensional criterion for matrix failure in compression becomes identical to his 1973 criterion. 

The coefficient of longitudinal influence, if, can be determined from shear tests with varying degrees 
of transverse compression. In the absence of biaxial test data, if can be estimated from the longitudinal 
and transverse shear strengths, as proposed by Puck: 





7 c eosa 0 (tan 2 a 0 sin a 0 + cos a 0 ) 


=> // 
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S L cos 2 a 0 
Y c cos 2 a 0 


( 11 ) 


Determination of the Angle of the Fracture Plane 

The angle of the fracture plane for a unidirectional laminate loaded in transverse compression is a 
material property that is easily obtained from experimental data. However, under combined loads, the 
angle of the fracture plane is unknown. The correct angle of the fracture plane is the one that maximizes 
the failure index, Ff in Eq. 4. In the present work, the fracture angle is obtained by searching for the 
maximum of the failure index (Eq. 4) within a loop over the range of possible fracture angles: 0 <a<a 0 . 
A graphical representation of matrix failure envelopes at various fracture angles is shown in Fig. 3 for a 
unidirectional E-Glass/LY556 composite subjected to transverse compression and in-plane shear. As seen 
in the figure, the fracture angle that maximizes the FI for small transverse stresses is cc= 0°. When the 
applied transverse stress a 22 has a magnitude equal to approximately 2/3 of the transverse compressive 
strength, the angle of the critical fracture plane switches from a= 0° to a=40°, and then rapidly 
increases to cc=a 0 , the angle of fracture for uniaxial transverse compression. 


T 12 , MPa 



Figure 3. Matrix failure envelopes for a typical unidirectional E-glass/epoxy lamina subjected to in- 
plane transverse compression and shear loading. 
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Criterion for Matrix Failure Under Transverse Tension (g^O) 


A failure criterion to predict matrix cracking under the presence of both in-plane shear and transverse 
tensile stresses should represent the “in-situ” effect occurring in laminated composites. The in-situ effect, 
originally detected in Parvizi’s 19 tensile tests of cross-ply glass fiber reinforced plastics, is characterized 
by higher transverse tensile and shear strengths of a ply when it is constrained by plies with different 
fiber orientations in a laminate, when compared with the strength of the same ply in a unidirectional 
laminate. The in-situ strength also depends on the number of plies clustered together, and on the fiber 
orientation of the constraining plies. 

The orientation of the constraining plies and the number of plies clustered together also affect the crack 
density and the stiffness reduction of the cracked ply. Wang’s 20 tests of [0/90 n /0] (n=l,2,3,4) 
carbon/epoxy laminates have shown higher crack densities for thinner 90° layers. The reduction of the 
elastic properties of a cracked ply is normally predicted using elastic analyses of cracked plies, 12, 21 or 
Continuum Damage Models. 22 ' 25 

The in-situ effect is illustrated in Figure 4, where the relation between the in-situ transverse tensile 
strength and the total thickness of the 90° plies clustered together is represented. 



Figure 4. Transverse tensile strength as a function of number of plies clustered together, with models 
from Dvorak 26 based on experimental data from Wang 20 . 

Accurate in-situ strengths are necessary for any stress-based failure criterion for matrix cracking in 

OH OO OQ 01 o OQ 

constrained plies. Both experimental ’ ’ and analytical methods ’ ’ have been proposed to 

determine the in-situ strengths. In the following sections, the in-situ strengths are calculated using 
fracture mechanics solutions for the propagation of cracks in a constrained ply. 


Fracture Mechanics Analysis of a Cracked Ply 

The failure criterion for predicting matrix cracking in a ply subjected to in-plane shear and transverse 
tension proposed here is based on the fracture mechanics analysis of a slit crack in a ply, as proposed by 
Dvorak and Laws. 26 The slit crack represents a manufacturing defect that is idealized as lying on the 1-3 


plane, as represented in Figure 5. It has a length 2a 0 across the thickness of a ply, t. Physically, this crack 
represents a distribution of matrix-fiber debonds that may be present in a ply as a consequence of 
manufacturing defects or from residual thermal stresses resulting from the different coefficients of 
thermal expansion of the fibers and of the matrix. Therefore, the slit crack is an “effective crack,” 
representing the macroscopic effect of matrix-fiber debonds that occur at the micromechanical level. 20 




Figure 5. Slit crack geometry (after Dvorak 26 ). 


Plane stress conditions are assumed. The transverse tensile stress o 22 is associated with mode I loading, 
whereas the in-plane shear stress Oi 2 is associated with mode II loading. The crack represented in Fig. 5 
can grow in the 1 (longitudinal, L) direction, in the 3 (transverse, T) direction, or in both directions. 

The components of the energy release rate for the crack geometry represented in Fig. 5 were 
determined by Dvorak and Laws. 26 For mixed-mode loading, the energy release rate for crack growth in 
the T and L directions, G(T) and G(L ), respectively, are given by: 


G(T) = ^ (if, A ° 22 <7 22 + Vn A 44 4 ) 
G(L) = ^ (f; A ° 22 (J 22 + A ° 44 4 ) 


( 12 ) 


where it can be observed that the energy release rate G(L) for longitudinal propagation is a function of 
the transverse slit size a 0 and that it is not a function of the slit length in the longitudinal direction a q . 

The parameters i=I f II in Eq. 12 are the stress intensity reduction coefficients for propagation in the 
transverse direction, and the parameters i=I,II are the reduction coefficients for propagation in the 
longitudinal direction. These coefficients account for the constraining effects of the adjoining layers on 
crack propagation: the coefficients are nearly equal to 1.0 when 2a 0 «t , and are less than 1.0 when a^t. 
Experimental results 28 have shown an increase in the in-situ transverse tensile strength of [±0/9O n ] s , 
0=0°, 30°, 60°, laminates for increasing stiffness of adjoining sublaminates ±0. This implies that the value 
of the parameter r/ x decreases with increasing stiffness of adjoining sublaminates. Considering that a 
transverse crack can promote delamination between the plies, Dvorak and Laws 26 suggested that the 
effective value of r/ x can be larger than obtained from the analysis of cracks terminating at the interface, 
and suggested the use of 7) x = £=\. 

The parameters A° jy are calculated as: 26 



V 


2 \ 


21 


"i y 


(13) 


9 


The mode I and mode II components of the energy release rate can be obtained for the T-direction 
using Eq. (12) with r/, = \ : 


Gl (T) = ^A\X 2 

(14) 

Gii(T) = ^A° 44 t?2 


The corresponding components of the fracture toughness are given as: 


G k (T)=^A\MJ 

(15) 

G,„(T) = 


where Y? and are the in-situ transverse tensile and shear strengths, respectively. 

For propagation in the longitudinal direction, the mode I and mode II components of the energy release 
rate are: 


G,(L) = ^ a » 22 <4 
G„(L) = ^ A° 44 4 

and the components of the fracture toughness are: 

G k (L) = ^A\Ml 

G, IC (L) = ^A”Mf 


(16) 


(17) 


Having obtained expressions for the components of the energy release rate and fracture toughness, a 
failure criterion can be applied to predict the propagation of the slit crack represented in Figure 5. Under 
the presence of both in-plane shear and transverse tension, the critical energy release rate G c depends on 
the combined effect of all microscopic energy absorbing mechanisms such as the creation of new fracture 
surface. Relying on microscopic examinations of the fracture surface, Hahn 30 observed that the fracture 
surface topography strongly depends on the type of loading. With increasing proportion of the stress 
intensity factor K Ih more hackles are observed in the matrix, thereby indicating more energy absorption 
associated with crack extension. Therefore, Hahn proposed a mixed mode criterion written as a first- 
order polynomial in the stress intensity factors Kj and K n . Written in terms of the mode I and mode II 
energy release rates, the Hahn criterion is 


0 - g\ 


'ML + g ML + M 0 =lf i=TtL 

G IC (i ) *G IC (i) G IIC (i) 


(18) 
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where the material constant g is defined identically from either Eq. 14 or 17 as: 


S = 


G IC A 0 , 


J lie 


22 


f y t ' 


A 0 , 


44 


V *is J 


( 19 ) 


A failure index for matrix tension can be expressed in terms of the ply stresses and in-situ strengths Y? 
and S l is by substituting either Eqs. 14-15 or 16-17 into the criterion in Eq. 18 to get: 
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( 20 ) 


The criterion presented in Eq. 20, with both linear and quadratic terms of the transverse normal stress 
and a quadratic term of the in-plane shear stress, is similar to the criteria proposed by Hahn 30 , Liu 31 (for 
transverse tension and in-plane shear), and Puck 5 . In addition, if g= 1 Eq. 18 reverts to the linear version 
of the criterion proposed by Wu and Reuter 32 for the propagation of delamination in laminated 
composites: 



J Ic 



= 1 


( 21 ) 


Furthermore, using g= 1, Eq. 20 reverts to the well-known Hashin criterion 9 for transverse matrix 
cracking under both in-plane shear and transverse tension, where the ply strengths are replaced by the in- 
situ strengths: 


fi m = 
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<1 


( 22 ) 


Application to thick embedded plies 

A thick ply is defined as one in which the length of the slit crack is much smaller than the ply 
thickness, 2 a 0 « t , as illustrated in Figure 6. The minimum thickness for a thick ply depends on the 
material used. For E-glass/epoxy and carbon/epoxy laminates, Dvorak and Laws 26 calculated the 
transition thickness between a thin and a thick ply to be approximately 0.7 mm, or about 5 to 6 plies. 




Figure 6. Geometry of slit crack in a thick embedded ply subjected to tension and shear loads. 
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For the geometry represented in Fig. 6, the crack can grow in the transverse or in the longitudinal 
direction. Comparing Eqs. 14 and 16, however, indicates that the energy release rate for the crack slit is 
twice as large in the transverse direction as it is in the longitudinal direction. Since Eq. 14 also indicates 
that the energy release rate is proportional to the crack length 2a 0 , the crack will grow unstably in the 

transverse direction. Once the crack reaches the constraining plies, it can propagate in the longitudinal 
direction, as well as induce a delamination. 

Crack propagation is predicted using equation (20), and the in-situ strengths can be calculated from the 
components of the fracture toughness as: 


(23) 


It can be observed from Eqs. 23 that the in-situ strengths of thick plies Y? and S L l s are functions of the 
toughnesses G Ic (T ) and G IIc (T ) of the material and the size of the material flaw size, 2a 0 . Therefore, the 
in-situ strengths for thick plies are independent of the ply thickness, as has been observed by Dvorak 26 
and Leguillon 33 , and as was shown in Fig. 4. 

Application to Thin Embedded Plies 

Thin plies are defined as having a thickness smaller than the typical defect, t < 2 a 0 , so the slit crack 
represented in Figure 5 extends across the entire thickness t of the ply, as represented in Figure 7. 



Figure 7. Geometry of slit crack in a thin embedded ply. 

In the case of thin plies, crack defects can only grow in the longitudinal (L) direction, or trigger a 
delamination between the plies. The in-situ strengths can be calculated from the components of the 
fracture toughness as: 


(24) 


where it can be observed that the in-situ strengths are inversely proportional to 4t . The toughnesses 
G Ic (L) and G IIc (L ) can be assumed to be the values measured by standard Fracture Mechanics tests, 
such as the DCB for mode I and the ENF test for mode II. Using Eqs. 23 and 24, Dvorak and Laws 26 
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obtained a good correlation between the predicted and experimentally obtained in-situ strengths of both 
thick and thin 90° plies in [0/90 n /0] laminates, as was shown in Fig. 4. 


Application to Unidirectional Laminates 

The fracture of a unidirectional specimen is taken as a particular case of a thick ply. 26 The defect size is 
2 a 0 , as for thick embedded plies. However, in the absence of constraining plies, the critical initial slit 
crack is located at the surface of the laminate. For tensile loading, the crack can be located at the edge of 
the laminate, which increases the energy release rate when compared with a central crack. In the case of 
shear loading, there is no free edge, so the crack is a central crack, as shown in Fig. 8. Crack propagation 
for unidirectional specimens subjected to tension or shear loads is obtained from the classic solution of 
the free edge crack. 34 



Figure 8. Unidirectional specimens under transverse tension and shear. 




where Y T and S L are the material tensile and shear strengths as measured from unidirectional laminate 

tests. Note that the in-situ strengths for thick plies can be obtained as a function of the strength of 
unidirectional laminates by substituting Eqs. 25 into 23: 


Yl = l.\l42Y T 

St = V2 S L 


The failure criterion for unidirectional plies under in-plane shear and transverse tension is represented 
in Eq. 20. The toughness ratio g for a thick laminate can also be calculated in terms of the unidirectional 
properties by substituting Eqs. 25 into Eq. 19 to get: 


g = 



A° 

1 . 12 2 2|1 


( y t V 


44 


V 


(27) 


LaRC03 Criteria for Fiber Failure 

Criterion for Fiber Tension Failure 

Based on the results of his experimental investigations, Puck 5 recommends using a fiber failure 
criterion that is based on the “effective strain” acting along the fibers. However, Puck also reports that 
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the difference between using the effective strain and the longitudinal stress is insignificant. 6 The LaRC03 
criterion for fiber tension failure is a non-interacting maximum allowable strain criterion that is simple to 
measure and is independent of fiber volume fraction and Young’s moduli. Consequently, the LaRC03 
failure index for fiber tensile failure is: 

LaRC03 #3 FI f = % < 1 (28) 


Criterion for Fiber Compression Failure 

Compressive failure of aligned fiber composites occurs from the collapse of the fibers as a result of 
shear kinking and damage of the supporting matrix (See for instance the representative works of Fleck 35 
and Soutis 36 , and Shultheisz’ 37 review of micromechanical compressive failure theories). Fiber kinking 
occurs as shear deformation leading to the formation of a kink band. 

Argon 38 was the first to analyze the kinking phenomenon. His analysis was based on the assumption of 
a local initial fiber misalignment. The fiber misalignment leads to shearing stresses between fibers that 
rotate the fibers, increasing the shearing stress and leading to instability. Since Argon’s work, the 
calculation of the critical kinking stress has been significantly improved with a more complete 
understanding of the geometry of the kink band as well as the incorporation of friction and material 
nonlinearity in the analysis models. " ’ 

Several authors 40, 41 have considered that misaligned fibers fail by the formation of a kink band when 
local matrix cracking occurs. Potter et al. 42 assumed that additional failure mechanisms that may occur 
under uniaxial longitudinal compression, such as crushing, brooming, and delaminations, are essentially 
triggered by matrix failure. 

In the present approach, the compressive strength X° is assumed to be a known material property that 
can be used in the LaRC03 matrix damage criterion (Eq. 4) to calculate the fiber misalignment angle that 
would cause matrix failure under uniaxial compression. 



Figure 9. Imperfection in fiber alignment idealized as local region of waviness. 


Calculation of Fiber Misalignment Angle 

The imperfection in fiber alignment is idealized as a local region of waviness, as shown in Fig. 9. The 
ply stresses in the misalignment coordinate frame m shown in Fig. 4 are 

<7™ = cos 2 <p <J n + sin 2 (p <j 22 + 2sin^cos^ r 12 
<J 22 = sin 2 cp<J n +cos 2 (p<J 22 -2s\n(pcos(pr n 

yZy) 

t™ 2 =-sin^cos^cr 11 + sin^cos^ cr 22 

+ (cos 2 (p - sin 2 (p ) t I2 
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At failure under axial compression, <J n = -X c , and <r 22 =T U = 0 . Substituting these values into Eq. 
29 gives <j ™2 = -sin 2 ^ c X c and r 12 = sinq) c cosq) c X c , where the angle qf is the total misalignment 
angle for the case of axial compression loading only. 

To calculate the failure index for fiber kinking, the stresses <T 22 and r 12 are substituted into the 
criterion for matrix compression failure (Eq. 4). The mode of failure in fiber kinking is dominated by the 
shear stress r 12 , rather than by the transverse stress a 2 2 • The angle of the fracture plane is then equal to 0°, 
and riff = 0 . The matrix failure criterion becomes 

T^r = X c (sin (p c cos qf - r/ L sin 2 (p c ) = Sf s (30) 


where is the in-situ longitudinal shear strength defined in Eq. 23 for thick plies and in Eq. 24 for thin 
plies. Solving for qf leads to the quadratic equation: 
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The smaller of the two roots of Eq. 31 is 
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(32) 


Note that if if were neglected and q> were assumed to be a small constant angle, Eq. 30 would give 


: ? which is the expression derived by Argon 38 to estimate the fiber misalignment angle. 




The total misalignment angle q> can be decomposed into an initial (constant) misalignment angle q)° that 
represents a manufacturing imperfection, and an additional rotational component qf that results from 
shear loading. The angles qf and qf can be calculated using small angle approximations and Eqs. 29 
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where the absolute value of the shear t 12 | is used because the sign of (p is assumed to be positive, 

regardless of the sign of the shear stress. Recalling that (p= (ft + (ft, it is now possible to solve Eqs. 33 
for (p in terms of <p c . 


<P = 


G 2 +(G n -X c )(p c 


G n 3" 


'22 


(34) 
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Fiber compression failure by formation of a kink band is predicted using the stresses from Eq. 29 and 
the failure criterion for matrix tension or matrix compression. For matrix compression (<r^ <0), the 
criterion is the Mohr-Coulomb criterion given in Eq. 4, with a=0° and = 0 . The criterion for fiber 
kinking becomes: 


LaRC03 #4 


FL 


*12 


+ 7J (J 


22 


st 


<1 


(35) 


For fiber compression with matrix tension, the transformed stresses of Eq. 29 are substituted into the 
matrix tensile failure criterion given in Eq. 20 to yield the following criterion for fiber kinking: 


LaRC03 #5 
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(36) 


It is interesting to note that for <j 22 = 0 , the fiber failure criterion in Eq. 35 becomes 


FI f 


-Cr il . p!2 1 

X c k s' 


<1 


with k — t - — > 1 

1 -<p 2 -2rj (p 


(37) 


The linear interaction between (7 n and r 12 in Eq. 37 is identical to the form used by Edge 7 for the WWFE. 
For T300/914C, Edge suggests using an empirical value of k= 1.5. However, Edge also indicates that 
other researchers have shown excellent correlation with experimental results with k= 1. Using the WWFE 
strength values of A 7 =900 MPa, .V /_ 80 MPa, Y c = 200 MPa, an assumed fracture angle in transverse 
compression of 53°, and Eqs. 11 and 32, we get: 7^=0.304, qf= 5.3°. With the approximation (p^(p C , Eq. 
37 gives ^=1 .07. 


LaRC03 Criterion for Matrix Damage in Biaxial Compression 

In the presence of high transverse compression combined with moderate fiber compression, matrix 
damage can occur without the formation of kink bands or damage to the fibers. This matrix damage mode 
is calculated using the stresses in the misaligned frame in the failure criterion in Eq. 4, which gives: 


LaRC03 #6 


fi m = 


f mT \ 2 
L eff 


+ 


f mL \ 2 
T eff 


<1 


(38) 


where the effective shear stresses and are defined as in Eq. 6, but in terms of the in-plane 
stresses in the misalignement frame, which are defined in Eqs. 29. 

= (-CJ 22 cos a (sma-r/ T cos a)) 

L / II, \ (39) 

T™ f f = ycos a (|tj 2 I + 7j a 21 cos ayj 


As for all matrix compressive failures herein, the stresses r ^ and are functions of the fracture 
angle a , which must be determined iteratively. 
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Verification Problems 


Example 1. Unidirectional 0° E-glass/MY750 epoxy 

All of the failure modes represented by the six LaRC03 criteria (Eqs. 4, 20, 28, 35, 36 and 38, and 
summarized in the Appendix) can be represented as a failure envelope in the On-o 2 2 plane. An example is 
shown in Fig. 10 for a 0° E-glass/MY750 epoxy lamina. Puck’s analysis results, which showed the best 
correlation with experimental results in the WWFE 3 , is also shown for comparison. In the figure, there is 
good agreement between LaRC03 and Puck in all quadrants except biaxial compression, where LaRC03 
predicts an increase of the axial compressive strength with increasing transverse compression. Testing for 
biaxial loads presents a number of complexities, and experimental results are rare. However, Waas et 
al. 43 report a number of references in which multiaxial compression was studied by superposing a 
hydrostatic pressure in addition to the compressive loading. For all materials considered, there is a 
significant increase in compressive strength with increasing pressure. In particular, the results of 
Wronsky and Parry 44 on glass/epoxy show a strength increase of 3.3 MPa per MPa of hydrostatic 
pressure, which gives an actual strength increase of 4.3 MPa per MPa of applied transverse biaxial stress. 
More recently, Sigley et al. 45 found 32% to 71% increases in compressive strength per 100 MPa 
superposed pressure. The results of Wronsky and Sigley can be compared qualitatively to the 4.3 
MPa/MPa increase in compressive strength for the plane stress failure envelope in Fig. 10. 


C 22 , MPa 



Figure 10. Biaxial Gn-G 2 2 failure envelope of 0° E-glass/MY750 epoxy lamina. 


Example 2. Unidirectional composite E-Glass/FY556 

A comparison of results from various failure criteria with the experimental results in the (J 22 -T 12 stress 
plane obtained from the WWFE 46 is shown in Fig. 1 1. It can be observed that within the positive range of 
< J 22 , all the quadratic failure criteria and LaRC03 give satisfactory results. Since the Maximum Stress 
criterion does not prescribe interactions between stress components, its failure envelope is rectangular. 
The most interesting behavior develops when o 2 2 becomes compressive. Hashin’s 1973 criterion gives an 
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elliptical envelope with diminishing r 12 as the absolute value of compressive a 2 2 increases, while the 
experimental data shows a definite trend of shear strength increase as a 22 goes into compression. 

The envelope for Hashin’s 1980 criteria was calculated using a transverse strength obtained from Eq. 
10, and it provides a modest improvement in accuracy compared to the 1973 criterion. Of the criteria 
shown in Fig. 11, Sun (Eq. 1), LaRC03 #1 (Eq. 4), and Puck (Eq. 2) capture the shear strength increase at 
the initial stage of compressive <j 22 . The failure envelope for Sun’s criterion (Eq. 1) was calculated using 
rf= 0.336 from Eq. 11. The results indicate a significant improvement over Hashin’s criteria. An even 
better fit would have been achieved using a higher value for rf. Puck’s envelope, which was extracted 
from the 2002 WWFE 6 , appears to be the most accurate, but it relies on fitting parameters based on the 
same test data. The LaRC03 curve uses the stiffnesses and strengths shown in Table 2, an assumed 
a 0 = 53°, and no other empirical or fitting parameter. In the case of matrix tension, Puck’s predicted 
failure envelope is nearly identical to LaRC03 #2. 

Table 2. Properties of E-Glass/LY556 46 ( MPa) 


E n 

E 22 

Gu 

V /2 

y t 

yC 

s L 

53,480 

17,700 

5,830 

0.278 

36 

138 

61 


WWFE test ‘02 Tl2 ’ MPa 



Figure 11. Failure envelopes and WWFE test data 2 for unidirectional composite E-Glass/LY556. 
Example 3. Cross-ply laminates 

Shuart 47 studied the compression failure of [±0\ s laminates and found that for #<15°, the dominant 
failure mode in these laminates is interlaminar shearing; for 15°<#<50°, it is in-plane matrix shearing; 
and for #>50°, it is matrix compression. Fiber scissoring due to matrix material nonlinearity caused the 
switch in failure mode from in-plane matrix shearing to matrix compression failure at larger lamination 
angles. The material properties shown in Table 3 were used for the analysis. The angle a 0 was assumed 
to be a typical 53°. For other lamination angles, the fracture angle was obtained by searching numerically 
for the angle that maximizes the failure criterion in Eq. 4. The in-situ shear strength did not need to be 
calculated since it is given by Shuart 47 as Sf s =95.1 MPa . 

The results in Fig. 12 indicate that the predicted strengths using LaRC03 criteria correlate well with the 
experimental results. The compressive strength predicted using Hashin 1973 criteria is also shown in Fig. 
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12. For #<20° , the Hashin criteria result in an overprediction of the failure load because the criterion 
does not account for the effect of inplane shear on fiber failure. For lamination angles near 70°, the use of 
the Hashin criteria result in an underprediction of the failure load because the criteria do not account for 
the increase in shear strength caused by transverse compression. All of these effects are represented by 
the LaRC03 criteria, which results in a good correlation between the calculated and experimental values. 

Table 3. Properties of AS4-3502 47 ( MPa) 


E 11 

e 22 

Gn 

V12 


Y° 

^ is 

127,600 

11,300 

6,00 

0.278 

1045 

244 

95.1 



Lamination angle, degrees 


Figure 12. Compressive strength as a function of ply orientation for [±0] s AS4/3502 laminates. 


Concluding Remarks 

The results of the recently concluded World Wide Failure Exercise indicate that the existing 
knowledge on failure mechanisms needs further development. Many of the existing failure models could 
not predict the experimental response within a tolerable limit. In fact, differences of up to an order of 
magnitude between the predicted and experimental values were not uncommon. In this paper, a new set 
of six physically-based failure criteria devoid of empirical variables is proposed. The criteria for matrix 
and fiber compression failure are based on a Mohr-Coulomb interaction of the stresses associated with 
the plane of fracture. Failure envelopes for unidirectional laminates in the a 'n-<J 22 and g 22 -t 12 stress 
planes were calculated, as well as the compressive strength of cross-ply laminates. The predicted failure 
envelopes indicate good correlation with experimental results. The proposed criteria, referred to as 
LaRC03, represent a significant improvement over the commonly used Hashin criteria, especially in the 
cases of matrix or fiber failure in compression. 
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Appendix. Summary of LaRC03 Failure Criteria 
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Required Unidirectional Material Properties: £„ £ 2 , G 12 , w 12 , A r , X c , Y t , Y c , S L , G Ic (L), G Sc (L ) 

Optional Properties: a 0 ,r/ L 


The in-situ strengths for a thin embedded ply are: 
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Toughness ratio g is obtained from fracture mechanics test data or from the unidirectional properties: 
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The effective shear stresses for matrix compression are 
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= ^cos a (|r 12 | + rj L a 22 cos or) ^ 


The transverse shear strength is S T = Y c cos a 0 ( sina 0 + CQSa ° j 


n 


-1 


tan 2 a n 


tan 2a 0 

and ~ - Sis CQS 2a ° or from test data. a 0 = 53 ° or from test data. 

Y c cos 2 


AQ ^ LUS tA0 

The stresses in the fiber misalignment frame are 

o™ = cos 2 ^< 7 n +sin 2 (p O' 22 + 2 sin cos (p T n 
0 22 = sin 2 (po n + cos 2 (p o 22 - 2 sin cp cos cp T u 
t \2 = - sin cp cos (po n + sin cp cos cp o 22 + (cos 2 cp- sin 2 (p) t 12 


where cp - 


z n + (G n - X c ) (p c 

G\2 + &\\ — 


and 


cp = tan 
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